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Abstract 

Recently, Csaki, Kaloper and Terning (CKT) suggested that axion-photon oscillation in the inter- 
galactic medium can explain the observed dimming of distant type la supernovae. This mechanism 
works only if the initial axion flux is much smaller than the initial photon flux. We study several 
mechanisms for such initial axion production. The direct axion production in the supernovae, the 
photon-axion oscillation in the magnetic field of supernovae and in the magnetic field of their host 
galaxies are addressed. We find it likely that the initial axion flux is very small and therefore does 
not pose a problem to the CKT mechanism. 
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I. INTRODUCTION 



Type la Supernovae (SNe la) are very interesting objects in observational cosmology [|TJ. 
Recent observations of distant (z ~ 1) SNe la have raised this interest enormously because 
the data indicate that they are fainter than would be expected for a decelerating universe 
|2|, 0]. The dimming seems to be independent of the wave length of the emitted light. The 
standard interpretation of this result is that our universe is accelerating at present as a 
result of a non- vanishing cosmo logical constant. Several alternatives to this idea have been 
suggested. One idea suggests that the light from distant SNe la can be partially absorbed 
by the dust which may be present on the light path ||. Ideas involving extra dimensions 
have also been suggested f|. The effect of cosmic evolution on distant SNe la was studied 

in I. 

Recently, Csaki, Kaloper and Terning (CKT) [[j] proposed another explanation of the 
dimming due to oscillation of photons into very light axions in the background magnetic 
field of the intergalactic plasma. The statistical significance of this result and possible 
generalizations of the mechanism are discussed in Ref. ||. 

The intergalactic plasma produces an energy dependent photon-axion oscillation and thus 



a frequency dependent dimming ||, [10], [TT[. This is in contrast to the achromatic dimming 
caused by a non- vanishing cosmological constant. While in the future the energy dependence 
may be the best tool to discriminate between these two mechanicsms, the present data are 



not sensitive enough |T0| , pJ]j . Based on this fact, it is important to further investigate the 
CKT mechanism. 

In all these works |7|, |8|, || [K| (El] it was implicitly assumed that only photons are emitted 



from the supernova (SN) while axions are produced only via oscillation in the intergalactic 
medium. If this assumption failed the dimming would be less significant. As an example 
consider the case of polarized photons, where the initial axion flux is equal to the initial 
photon flux (with polarization state parallel to the external magnetic field). Then, since 
the photon to axion and axion to photon oscillation probabilities are equal (P 7 _> a = Pa^-y) 
the oscillation will cause no effect at all. More generally, in the case of unpolarized photons 
and/or randomly distributed background magnetic field, the photon intensity on Earth is 
given by 

J 7 (y) = 1,(0) - P^ a {y) [1,(0) - 2J o (0)] (1) 



2 



where y is the distance traveled by the photons, I 7 (0) [J o (0)] is the initial photon [axion] 
flux, I 1 {y) is the final photon flux and P 1 ^ a {y) is the oscillation probability. Clearly, if the 
initial axion flux is significant it must be taken into account. 

In this paper we estimate the initial axion flux for the relevant energy range. We consider 
axion production in the SN, axion conversion in the SN magnetic field and axion conversion 
in the magnetic field of the host galaxy. We find that it is unlikely that the initial axion 
flux is significant. We thus conclude that neglecting the initial axion flux does not pose a 
problem to the axion-photon mixing explanation of the dimming. 



II. PHOTON- AXION MIXING FORMALISM 



In a background magnetic field axions and photons oscillate into each other [T2|, [L3[ due 
to the axion-photon coupling 

tint = j^E ■ B, (2) 

where a is the axion field, E is the electric field of the photon, B is the background magnetic 
field and the mass scale M characterizes the strength of the axion-photon interactions. In 
Ref. [[7] it was found that M = 4 x 10 11 GeV gives the best fit to the data and we will use 
this value in the analysis. In order to have mixing a transverse external magnetic field is 
needed. Moreover, only the plane wave photon polarization state parallel to the external 
field can mix with the axion. 

Photon-axion oscillation is induced by the following mass-squared matrix 

\ — it /j, m ) 

where £ is the photon energy, 

" = (4) 

B ~ |B| is the projection of the magnetic field on the direction of the photon polarization 
and m is the axion mass. In Ref. 0, the authors found that m ~ lCT 16 eV is needed. For 
the purposes of this paper, as long as m 2 max(£ fi, u 2 ), we do not need the exact value of 
the axion mass (m ~ lCT 16 eV is small enough). The plasma frequency is given by 

u 2 p = e - , 5 
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where m e is the electron mass and n e is the number density of electrons in the plasma. We 
omitted the Euler-Heisenberg term |b| uj 2 eh = (7a/907r)(5 /B cr ) 2 £ 2 (where B cr = 4.41 x 
10 13 G is the critical field strength which contribute to since the magnetic fields 

we are considering are much smaller than B cr and the effect of this term is negligible. 
For a constant magnetic field the oscillation probability of photons into axions is given 

by 

fi 2 . 2 (kx\ {fix) 2 ( sm(kx/2) \ 2 

where x is the distance traveled by the photon, k = 2n / ' Lq sc and the oscillation length is 
given by 

T (7) 

L Osc — I 9 = ■ { I ) 

tJ{uj 2 - m 2 f + Afi 2 S 2 
Note that the oscillation probability is bounded by 

Pl <^f. (8) 

Of particular interest to us is the case where the constant magnetic field domain size (Ld om ) 
is small, namely when kLd om "C 1. Expanding in this small parameter one finds that the 
bound of Eq. (§) (with x = L dom ) is saturated, namely Pi = (/iL rfom ) 2 /4. In particular, 
in this case the conversion is independent of the axion mass, the photon energy and the 
plasma frequency. This is the reason that the CKT mechanism does not produce energy 
dependent dimming as long as the plasma frequency is small enough such that the above 



approximation is valid 1C 



Next we consider the case where the magnetic field is not constant. Within a set of sim- 
plifying assumptions (see the Appendix for the derivation) the total photon-axion conversion 
over many domains is given by 

P^ a (y) = l(l- e~y' L ^y) , (9) 
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such that y = NLd om is the total length traveled by the photons in iV domains and 

J- 1 decay j 

where P\ is given in (H). Using (H) we find the following bound 



(10) 



As discussed in the appendix there are several assumptions that enter the derivation of @. 
Still, we do not expect significant deviation from the bound ([11]) in more realistic cases. 

We also consider cases where the mean free path of photons, L m f p , is not very large. In 
such cases it takes an average of N = 3i? 2 /L^j p random steps to escape from a region of 
radius R. Assuming that the magnetic field is constant over a region of order L m f p , the 
conversion probability is given by Eq. ([]) with 



V — NL m f p — , L decay — (12) 

Using (JTTJ) we get the bound 

p^ a (y = R) < ^f^ ■ (is) 

Note that the bound does not depend on the photon mean free path. 

In the following we will make use of the bounds @, ( 11 ) and (pT3[). Of course, they are 



useful only when the characteristic scale set by the magnetic field (l//i) is much larger than 
the relevant scales in the problem (x, y and R). As it turns out this is indeed the situation 
in the cases we study. 

III. PHOTON AXION MIXING INSIDE SUPERNOVA IA 

We now turn to study axion-photon mixing inside the supernova and determine under 
which physical conditions there is an appreciable photon-axion oscillation probability. SNe 
la are identified by the variations of their apparent luminosity with time. The intrinsic 
peak luminosity of these SNe varies within 20 %, and based on the time variation even this 
variation can be corrected for. This is the reason that they are used as "standard candles" for 
various cosmological tests. The spectral analysis of SNe la shows the absence of hydrogen. 
Supernovae of this type are believed to result from the thermonuclear explosion of White 
Dwarfs (WDs) when their masses are very close to the Chandrasekhar limit |15, [16 . 



In order to estimate the axion-photon oscillation probability inside SNe la we consider 



a simple model W<n - The supernova is assumed to be a sphere of uniform density with an 
initial radius Rq ~ 10 9 cm that expands with an outer velocity of v = c/30 = 10 9 cm s _1 . 
Photons are emitted uniformly throughout the volume of the supernova and the energy 



distribution follows Planck's law [17, 18]. The peak luminosity occurs within 10 to 20 days 



after the explosion. The progenitor stars are magnetized objects. Some WDs are known to 



host magnetic fields of intensity ranging between 10 5 and 10 9 G. The magnetic fields during 
the supernova event are not very well understood. It is possible that due to the turbulent 
motions prior to the explosion the WD goes through a stage which might dramatically 
change the intensity of its initial magnetic field. In particular, if the effect is to increase the 
magnetic field its value can be as large as B ~ 10 11 G [|19|| . After the explosion of the WD, 
the ejecta undergo a large expansion. In a homologous expansion all components of the field 
decrease like 

B_ = {R^_ ^ (Rol 

B (R) 2 ~ {vtf ' 1 1 

where the approximation is valid for vt 3> Rq. This background magnetic field can be 

chaotic which is a likely result of the turbulent motion prior to the explosion. 

We first consider the production of axions in the SN envelope due to photon conversion 

in its magnetic field. The photon mean free path at peak luminosity is not well known. 

Ref. | 18fl , for example, estimates that L m f p is in the range 10 6 — 10 14 cm. The SN radius 



at that time, which is of order 10 15 cm, is larger than the mean free path. The conversion 
probability at peak luminosity is therefore bounded by Eq. (f[3|). Using ([14]), we get 

P _ < ? (mv „ 3 x 10 - ( * y fi^s-y (15) 

1 a -8\MvtJ V10 n G/ \ t J K ' 

where for the last step we used M = 4 x 10 11 GeV, Rq = 10 9 cm and v = c/30. We learn 
that for values of t > 10 days and B Q < 10 11 G the conversion probability is negligible. Our 
conclusions would not hold in the unlikely case that initial magnetic fields in SNe la reach 
values significantly higher than 10 14 G. 

Let us now discuss other processes by which axions could in principle be produced inside 
the SN. It is well known that axion production is potentially large in the first stage of type II 
supernovae explosions [EQ]. One would then like to check whether a similar situation occurs 
in the type la SNe explosions. However, this is not the case since the temperature and the 
density of matter in the plasma of SNe la at its peak luminosity are relatively very small. 
Let us elaborate on this. 

The potentially relevant processes are the ones that occur at the typical energy scale 
of a few eV. The axion luminosity due to these processes scales as n p B n^T s , where n# is 
the baryon density, n e the electron density, T the temperature, and p, r, and s are process 
dependent nonnegative numbers. We consider the following processes: axiorecombination 
- i.e. the reaction, where axions are emitted when electrons and ions form a bound state; 



axion bremsstrahlung from the electrons and nuclei in the plasma; Compton scattering 
with axion instead of a photon in the final state; and the axion emission from electrons 
via the conversion of longitudinal plasmons. The luminosity due to the axiorecombination 
scales as r2,£?7, e T 3 / 2 |JI|], the luminosity due to the bremsstrahlung from electrons scales as 



nsn e T 5//2 |21J • The axion bremsstrahlung emission from nuclei in the eV energy range is 
suppressed by additional powers of m e jra^. The luminosity due to the Compton scattering 
scales as n e T 6 |20|, and the luminosity due to the plasmon conversion scales as n^r||22|. The 



density of plasma in SNe la at the time of peak luminosity is of order 10 13 gem 3 and 
the temperature is about 2 x 10 4 K [TJ|]. To get a rough estimate of the axion production 



in the SNe la we compare these values with the typical density p ~ 10 13 gem 3 and the 
temperature of several MeV inside the core of SNe II [p0| |. Multiplying by the volume of the 



SNe explosion remnant and taking n e ^P? roughly constant we arrive at a suppression factor 
of about 10~ 40 — 10~ 50 . This expectation is borne out by a more detailed calculation, where 
it is found that the axion flux due to these processes in the SNe la remnant is suppressed 
by about 50 orders of magnitude compared to the photon flux. We conclude that axion 
production by these mechanisms can be safely neglected. 

IV. PHOTON AXION MIXING IN THE HOST GALAXY 

Next we consider photon axion conversion in the host galaxy. The details of the back- 
ground magnetic field of the galaxy are not well understood. We first assume that far away 
galaxies are similar to the better known nearby galaxies. The size of the galactic disk, 
which is also the coherence length of the smooth component of the magnetic field is of order 
1 kpc ~ 10 26 eV" 1 . For a typical galaxy, Ref. [^3j estimates a smooth toroidal component 
of the magnetic field of magnitude about 2 fiG and a random component with a magnitude 
of about 5 fiG and a coherence length of order 50 pc. The mean electron density in the host 
galaxy is taken to be n e = 0.03 cm -3 p3| , With these parameters we find that the mean 
free path of photons is much larger than 1 kpc and thus they basically do not interact after 
they leave the SN. For the smooth toroidal component of magnetic field, using B = 2 /iG, 
we have /i ~ 10~ 28 eV. Using (||) with x = Ld om = 1 kpc we conclude 

P 7 ^ a <10" 4 . (16) 
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The random component of the magnetic field is somewhat larger in strength. However, due 
to smaller coherence length its contribution to the conversion probability is likely to be even 
smaller. For example, using (|TT| ) with L dom = 50 pc, y = 1 kpc and B = 5 /iG we find 

P^ a < lO" 5 . (17) 

Far away host galaxies may have considerably different magnetic field than what we find 
in nearby galaxies. One or two order of magnitude enhancement of the magnetic field could 
result in a significant conversion. (In order for such conversion to be realized also the plasma 
frequency has to be small.) Yet, there is no well motivated reason to believe that this is the 
case. Thus it is likely that the conversion probability in the far away host galaxy is small. 

V. CONCLUSIONS 

We have studied axion production in type la supernovae and photon-axion oscillations 
in their magnetic fields and in the magnetic fields of their host galaxies. If such axion 
production is significant it can disfavor the photon-axion oscillation explanation of the 
observed dimming of light from far away type la supernovae since then the dimming is 
less significant, and it may also become energy dependent. We found that all these effects 
are very small. Our result is robust in the sense that it does not rely on many unknown 
parameters of the SNe and their host galaxies. For example, we did not have to use exact 
values of the plasma frequency or the photon mean free path. We had to use only the rough 
size of the star, its temperature and its magnetic field. We therefore conclude that axion 
production in the supernovae and their host galaxies has a negligible effect on the CKT 
explanation of the supernovae dimming. 
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APPENDIX A: CONVERSION IN VARYING BACKGROUND 



In this appendix we discuss photon to axion conversions in the varying magnetic field 
background. To simplify the discussion it is assumed that photons and axions traverse N 
domains of equal length. In each domain the magnetic field B is assumed to be homogeneous 
with a discrete change in B from one domain to the other, i.e. the situation is far from 
adiabatic conditions. The component of the magnetic field perpendicular to the direction 
of flight is assumed to have random orientation but is equal in size in each domain. The 
domains can be (not necessarily) separated by an interval where the conversion is negligible. 

Consider an initial state Ci(0)|7i) + c 2 (0)|7 2 ) + c a (0)\a). The two photon polarization 
states, 171,2}) correspond respectively to photons parallel and perpendicular to the magnetic 
field in the first domain. The initial photon and axion fluxes are then I 7 (0) ~ | ci (0) | 2 + 
|c 2 (0)| 2 and I a (0) ~ |c a (0)| 2 respectively. In the n-th domain the magnetic field is tilted by 
an angle 9 n compared to the magnetic field in the first domain 

|7f) =c n | 7l )+s n |7 2 ), = -s n \ lx ) + c n | 72 ) , (Al) 

or 

c x {y) = c n c\(y) - s n cl(y) , c 2 (y) = s n c\(y) + c n cl(y) , (A2) 

where c n = cos# n and s n = sm8 n . The magnetic field mixes photons polarized parallel to 
the magnetic field with the axions according to the mass matrix given in (^). We denote 
the transition elements for one domain t 7 ^ 7 = (7|f(£dom)|7|f(0)) and similarly for t^ a . 
The values of the transition elements are equal in each domain since the magnitude of the 
magnetic field has been assumed to be the same everywhere. The transition probability for 
photon to axion oscillation in one domain is Pi = |t 7 ^ a | 2 with P\ given in (|5p, while the 
survival probability of photon is 1 — P\ = |t 7 ^ 7 | 2 . At the end of the n-th domain the photon 
and axion fluxes are 

J 7 (n + 1) ~ (1 - P lC 2 n )\ Cl (y n )\ 2 + (1 - P lS 2 n )\c 2 (y n )\ 2 + Pi|c a (i/ n )| 2 + . . . (A3) 
I a (n + 1) ~ Pic 2 |ci(y n )| 2 + Pi S 2 |c 2 (i/ n )| 2 + (l-Pi)|c a (t/ n )| 2 + ... (A4) 

where dots represent terms that are proportional to cos# n , sm9 n or cos0 n sin# n . We have 
defined y n = (n — l)Ld om . Namely, the coefficients Ci j2ja are taken at the beginning of the 
n-th domain. 
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Next we assume that the transition probability in one domain is small i.e., Pi < 1, 
and that the direction of the magnetic field is random, i.e., 8 n is a random variable with 
$n+i — ®n ~ In this limit c 2 and s 2 can be replaced by their average value 1/2, while 

c„, s n and c n s n are averaged to zero. Using the fact that I 7 (n) ~ |ci(2/ n )| 2 + |c 2 (|/ ri )| 2 and 
I a ~ l c a(?/n)| 2 we arrive at 

1 - \P X Pi 
\Px I- Pi 
1 



/ 7 (n+l) 
Mn + l) 



i" 7 (n) 
/„fn) 



(A5) 



3 p 
3 p 



n+1 
n+1 



1+2 
3^3 



3 p 



3 p 



n+1 
n+1 



•/ 7 (o) 
.4(0) 



We can then use the fact that the number of domains is large and replace (1 — 3Pi/2) ra+1 
with the limiting expression exp [— (3Piy) / (2Ld om )} to arrive at the final expressions 



L f {y) = A(o) - P 7 ^(/ 7 (o) - 2/ a (o)) , 

I a (y) = / a (0)+P 7 _ a (/ 7 (0)-2J o (0)), 



(A6) 
(A7) 



with 



I- 



-y/Ldec 



' CaV ) > ( A8 ) 

and L decay = (2L dom )/(3Pi). In cases where the oscillation length is much larger than the 
domain size the conversion probability Pi saturates the upper bound Pi = (/iLd om ) 2 /4 and 
we arrive at the expression Ldecay = 8/(3yU 2 Ld om ) given also in |?J. 
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